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The symmetric collective states of an atomic spin ensemble (i.e., many-body states that are invari- 
ant under particle exchange) are not preserved by decoherence that acts identically but individually 
on members of the ensemble. We develop a class of collective states in an ensemble of A'^ spin-1/2 
particles that is invariant under symmetric local decoherence and find that the dimension of the 
Hilbert space spanned by these collective states scales only as A'^'^. We then investigate the open sys- 
tem dynamics of experimentally relevant non-classical collective atomic states, including Schrodinger 
cat and spin squeezed states, subject to various symmetric but local decoherence models. 
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Many fields, including precision metrology [T, ^5], quan- 
tum information science [Sj i4j (5n6j and quantum optical 
models of condensed matter phenomena [II [H [9], uti- 
lize collective spin degrees of freedom in an ensemble of 
atoms or ions. An attractive feature of such ensembles is 
that one can achieve exquisite control over their quantum 
dynamics, even in practice A second benefit is 

that interesting many-body effects can be realized within 
a small, theoretically manageable, sub-Hilbert space of 
the entire spin ensemble's state space. Even though the 
joint Hilbert space for N spin-j particles is immense, 
dim(J^) — (2j'-fl)^, it is not uncommon to find dynam- 
ical symmetries that reduce the effective dimension of the 
ensemble without sacrificing the rich structure character- 
istic of many-body interactions and entanglement. 

Most previous work on dynamical symmetries in spin 
ensembles has involved the symmetric collective states 
[TU]: the sub-Hilbert space Ms C spanned by A^- 
body states that are invariant under the permutation 
of any two particle labels HijlV') — IV') i IV') G 
The dimension of this so-called "symmetric space" grows 
only linearly in the number of spin particles dim(^^) = 
{2Nj + l) <C (2j-t-l)^, yet exhibits phenomena of consid- 
erable interest, such as spin-squeezing [TT1|T2] and zero- 
temperature phase transitions [T . For to provide an 
accurate description of the ensemble's state, the dynam- 
ics must be generated by completely symmetric collective 
processes — symmetric indicating that the processes are 
permutation invariant, and collective meaning that they 
can be expressed in terms of total angular momentum 
operators. Symmetric states are accessed in practice by 
experiments that avoid particle addressability, such as 
by homogeneous coupling to the external fields used to 
control and observe the ensemble [5J [10] . 

Unfortunately many of the decoherence models most 
appropriate for large spin ensembles cannot be described 
as collective processes even when the decoherence acts 
identically on each particle. The {N + l)-dimensional 
symmetric Hilbert space Ms is not preserved under such 
decoherence, greatly complicating dynamical simulations 
of the ensemble. It has thus become common practice 



either to: (1) restrict the dynamics to very short-times 
when decoherence can be approximately neglected; or (2) 
model the decoherence as a collective process, even when 
doing so is not necessarily physically justified. In atomic 
spin ensembles, for example, a typical source of deco- 
herence comes from spontaneous emission, yet collective 
radiative processes only occur under specific conditions 
such as superradiance from highly confined atoms [ 13] 
and some cavity-QED or spin-grating settings [^. In 
many ongoing experiments, fiuorescence imaging can, to 
varying degree, distinguish between particles. 

In this Letter, we generalize the collective states of an 
ensemble of spin-1/2 particles (qubits) to include states 
that are preserved under symmetric — but not necessar- 
ily collective — transformations. Specifically, we gener- 
alize from the strict condition of complete permutation 
invariance to the broader class of states that are indistin- 
guishable across degenerate irreducible representations 
(irreps) of the rotation group. The representation the- 
ory of the rotation group has, of course, been utilized 
in a wide variety of contexts, such as to protect quan- 
tum information from decoherence by encoding it into 
degenerate irreps with the same total angular momen- 
tum [131 [TS]. Here, we utilize relevant aspects of the 
representation theory to obtain a reduced-dimensional 
description of quantum maps that act locally but iden- 
tically on every member of an ensemble of qubits. The 
dimension of the Hilbert space Mq spanned by our gener- 
alized collective states scales favorably, dim(jJ^) ~ N"^ , 
making it possible to address decoherence models rele- 
vant to a greater variety of experimental settings. To 
demonstrate the viability of our approach, we investi- 
gate several decoherence models applied to Schrodinger 
cat and spin-squeezed states of an atomic ensemble. 

Representations of the Rotation Group. — We consider 
an ensemble of N identical spin-1/2 particles, each de- 
scribed by the single-particle angular momentum oper- 
ators tr^") = ((Ti"\ (Ty"\ ai"^). The joint Hilbert space 
for the entire spin ensemble M — M'^^^ ® ■ ■ ■ (i) Jf'^^^ 
has dimension dim{M) — 2^, and arbitrary pure states 
of the ensemble can be expressed in the tensor prod- 
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|mi, 7712, . . . , TOjv) where 
D • • • ® I i , m^v) jv satisfy 



uct basis IV") = X^r 
the |mi, . . . , mjv) = 

ji"^|mi, . . . ,TOAr) = m„|mi, . . .,mN)- 

Each particle in the ensemble transforms separately 
under a rotation such that jV'') = (i?)]"^^]^)^ where 
^3(i?) is the spin-1/2 rotation operator parameterized 
by the Euler angles R — (a,/?, 7). Expressed in the 
tensor-product basis, the [^2 (i?)]®^ provide a reducible 
representation for the rotation group but can be decom- 
posed into irreducible components (irreps) 



(1) 



i=i 



via the total spin eigenstates j^|J, Af, i) = J(J 
l)\J,M,i) and X\J,M,i) = M|J,M,i 



with Jq = 



Y- For each to- 



|E„=i'^r and J = mod(f ,2), 
tal angular momentum J, the quantum number i 
1, . . . , dj^ distinguishes between the 
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N\{2J+1) 



{N/2 - J)\{N /2 + J +1)1 



(2) 



degenerate irreps with total angular momentum J 

Collective States. — In the "irrep basis," arbitrary pure 
states of the spin ensemble are expressed as \tp{t)) = 
M i'^JM,i{'t)\J,MTi), which still requires 2^ coeffi- 
cients. Of course, simply transforming to the irrep basis 
does not change the effective dimension of the Hilbert 
space, but it suggests a symmetry that can be exploited 
to do so. We generalize the definition of collective states 
to the sub-Hilbert space spanned by iV-particle states 
\'4'c) that are indistinguishable across the df^ degenerate 
irreps for each total angular momentum J. This intro- 
duces the symmetry cj^M,i = cj^M,i' , making it im- 
necessary to distinguish \J,M,i) over irrep label. Defin- 
ing effective kets |J, Af) on each total- J irrep, we then 
have IV'c) = m'^J,m\'^t^) ^^^^ the coefficients 



Cj,M{t) 




(3) 



Under this symmetry, dim(J«?) = {N + 2)^/4 (for N 
even) scales only quadratically with the number of parti- 
cles and therefore admits efficient representation. When 



describing decoherence and other open-system dynam- 
ics, such as continuous measurement, one may require 
the collective state density operator, which we define as 
PC = ®jPj to give 

Pcit) = E E PJ,M:J,M'itMM){J,M'\. (4) 

J M,M' 

These collective density operators restrict against coher- 
ences between irreps [17] , but still subsume the symmet- 
ric space previously considered for large spin ensembles 
|10j . Jifg C .y^c is spanned by the maximal angular mo- 
mentum manifold: cj^m = for J 7^ N/2. All \ips) G J^s 
can be obtained from the permutation-invariant state 
\N/2,N/2) <-> |i . . . , i) by a simple rotation ^^/^(i?) 



and are collective states since d 



N/2 
N 
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Collective State-Preserving Processes. — We next dis- 
cuss quantum maps that preserve collective states C : 
Jfc '^G- The super-operator 



PJ,M\J,M' f'llM 



(5) 



J M,M' 



will be closed on if its action //^^,^, = C\J, M) (J, M'\ 
can be expressed in the irrep basis in such a way that it 
does not distinguish between degenerate irreps. Opera- 
tors that transform simply with respect to the rotation 
group = E J E MM' ^M,M' I J' ^) ^' I ' including the 
collective angular momentum operators Ji and all collec- 
tive operators S — E^=i ^^"^ formed from s^"^ G su(2), 
satisfy this requirement directly. 

Unfortunately, not all dynamics can be expressed en- 
tirely in terms of collective operators; they might only be 
symmetric over local single-particle super-operators. 



C[s]p - 



N 
n=l 



(6) 



where /:(")[s(")]/5 = s(")/5(s("))^ for e su(2). Our 
main result is that any symmetric local map of the form 
Eq. ^ can be brought into the form of Eq. ^ and 
therefore preserves collective states. For the 5u(2) op- 
erators s = s ■ <T and = t ■ a, f^fj^fi can be con- 
structed as /mjv/' — ^' s{J:M,M') ■ t from the tensor 
9qr = E^=i4"V,^)('>^,^'l('^r"V whosc elements 
find to take the form [18] 
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FIG. 1: (color online) Decoherence of an initial Schrodinger cat state: (al-a2) time-dependent fidelity with the initial state for 
both symmetric local /^[o"-] and collective /^[J-] decoherence for different numbers of particles; (bl-b2) similar comparison for 
£[(Tz] versus £[jz]; (c) time-dependent populations of different total-J irreps for -C[i5"-]. 



where the indices q,r e { — ,-|-,2;,0} denote representing 
single-particle angular momentum operators in the basis 
<T+, (Tz, 1}. We have also defined the following sym- 
bols: M+ ^ M + I, M_ = M - 1 and = M, the 

reduced degeneracies a'}^ = J2j'=,j'^ 
coefficients 



j(r , and the tensor 



+ ^(JzpM)(J±M + 1) 
M 

±^{JtM){JtM-\) 



B 
D 



J,M _ 



J,M 



= ^{J + M){J-M) 

TViJ±M + 1){J±M + 2) 
v/(J + M+l)(J-Af + l) . 



(8) 
(9) 
(10) 

(11) 
(12) 

(13) 



The three terms in Eq. (|7| arise from two types of pro- 
cesses: (Term 1) transitions that occur between M levels 
within a single J irrep; and (Terms 2-3) transitions that 
couple neighboring irreps with A J = ±1. For brevity, we 
have not listed terms in Eq. ([t]) for q = (the identity 
operator); these terms as well as those listed here can be 
determined following the methods in Ref. |T5]. 

Reducing the effective dimension of the ensemble 
Hilbert space to 0{N^) should be adequate to allow sim- 
ulations of N > 100 on a desktop computer, and much 
more on specialized hardware. It is, however, worth men- 
tioning that the structure of the mapping in Eq. ^ sug- 
gests an opportunity for further reducing the effective 
dimension of the ensemble's Hilbert space. Since the dy- 
namics only couple neighboring irreps, if the initial state 
is confined to only a single irrep (for example a symmet- 
ric state) then it may well be possible to truncate distant 
irreps from the decomposition without introducing ap- 
preciable error into the ensemble dynamics. 



Modeling Decoherence. — As discussed in the introduc- 
tion, realistic decoherence models for an ensemble of spin 
particles are often described most aptly by a symmetric 
sum over local channels. Consider, for example, the open 
system dynamics governed by the master equation 



dm 

dt 



^~i[H,p{t)] + TC[s\p{t), 



(14) 



where H (and any measurements performed) are de- 
scribed by collective operators, but the decoherence in- 
volves the symmetric Linblad superoperator 

N 

C[s\p = E - 2 + p[st5](")) (15) 



of the form in Eq. As this decoherence model does 
not preserve symmetric states, it has been common prac- 
tice to consider instead the associated collective process 



SpS^ - ^ {s^Sj^ 



bs^s 



(16) 



with 5* = ^„s("'. 

To illustrate the difference between symmetric and col- 
lective decoherence models, we considered the open sys- 
tem dynamics of two representative problems. First, 
we compared the dynamics generated by the symmetric- 
local C[s\ versus collective C[S] Linblad master equa- 
tions applied to an initial Schrodinger cat state |'0(O)) = 
(If ,+f ) + If '-f )) /v^- Figure |l[a-b) depicts the fi- 
delity T{t) = (V'(O)|p(t)|'0(O) evolwd under Eq. ([l4l 
(with H ^ 0) for two different types of decoherence 
channels: Fig. [ljal-a2) compares the collective versus 
symmetric master equations with s = dz for A'^ = 10 
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FIG. 2: (color online) Time-evolution of the squeezing pa- 
rameter for a spin ensemble driven \yy H = — J^) 
subject to -C[(5"-] (dotted lines) and /I[J-] (dashed lines) with 
relative decoherence rates V — A/5,A, 5A. For comparison, 
the solid line denotes decoherence-free squeezing. 



population in each irrep block then undergoes its own 
squeezing, evidently making the dynamics more resistant 
to symmetric local decoherence than collective processes. 

Conclusion. — We have presented an exact reduced- 
dimensional description of quantum maps that act locally 
but identically on each member of an ensemble of qubits. 
Our results confirm that the common practice of model- 
ing decoherence in large atomic ensembles via a collective 
process (when not physically motivated) is unlikely to al- 
low for quantitative comparison between theory and ex- 
periment. The techniques developed here are expected to 
help achieve such agreement, and thus play an important 
role in quantum measurement, control and information 
theoretic settings where high-fidelity modeling is essen- 
tial to interpreting experimental data. This work was in- 
spired by discussions with Ivan Deutsch and Poul Jessen 
and was supported by the NSF (PHY-0652877) and the 
DOE under contract with Sandia National Laboratory 
through the NINE program (707291). 



and N = 100 particles, respectively; and Fig. [TJbl-b2) 
makes a similar comparison for s = a^- The examples 
we considered (including some not reported here) sug- 
gest symmetric local decoherence models can generate 
dynamics that are appreciably different from their col- 
lective analogs. This is perhaps not too surprising: for 
an initially symmetric state, collective decoherence mod- 
els C[S] confine the dynamics to only maximum- J irrep; 
symmetric local models C[s] do not necessary preserve 
the irrep decomposition of the initial state. Fig. [ijc) de- 
picts the norm of each total- J irrep block of the density 
operator Nj — tr[Pjp{t)] as a function of time for £[(T-] 
(Pj = I]mI^'^0(^'^^I)- The observation that small- J 
irreps are only minimally populated suggests that further 
model reduction by truncating the Hilbert space to only 
the largest J blocks could be beneficial. 

As a second example, we compared symmetric- 
local versus collective decoherence models applied to 
dynamically-generated spin squeezing under the counter- 
twisting Hamiltonian H = — iA( — J^) |TTj. We per- 
formed simulations by time-evolving Eq. ( 14 ) from the 
initial spin-coherent state |^,^) for N = 100 with 
C[a-] and £[J_]. Figure § depicts the time-dependent 
squeezing parameter = N{AJy)/{Jz)'^, each for F = 
A/5,A, 5A. Under the conditions we considered, sym- 
metric local decoherence wave evidently less destruc- 
tive to the squeezing dynamics than collective models. 
As observed for the cat-state dynamics, to large extent 
the main effect of symmetric-local decoherence is leak- 
age from the maximum J irrep. But since the driving 
Hamiltonian H involves only collective spin operators, 
the coherent dynamics decouple for different total J: the 
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